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ABSTRACT 

An imoroved Objective Analysis technique (Gendin, 1963) ia naeu along with an 
altimeter signal atatiatical model, an altimeter noise statistical model, an orbital 
model, and synoptic anriace current nap % in the POLYMODE-SDE area, to evaluate the 
performance of varioua observational strategies in catching the aesoscale variability 
at mid-latitndea . In particular, simulated repetitive nominal orbits of ERS-1, 

TOP EX, and SPOT/ POSEIDON are examined, -ur statistical models a*: consistent with 
previona in-»itu and remote-sensed results (Fu, 1983). Our own results show the 
critical i^ortance of the exists"**! of a subcycle, scanning in either direct; ^. 
Moreover, loug repeat cycles 020 ways) and short cross-track distant- (<300 kei) 
seem preferable, since they match aMsos'ale statistics. Another goal of onr study ia 
to prepare and discuss Sea-Surface Height assimilation in quaaigeostrophAc '-.dels. 

Our restored SSH maps are shown to meet that purpose, if /. n efficient extrapolation 
method or deep in-aitu data (floats) are used on the vortical to start and update the 
model . 

1. INTRODUCTION 

Active and passive satellite-flown sensors have been extensively slcv« to open 
new perspectives in geophysics.’ sciences; in particular, the altimeters of GEOS 3 and 
SEASAT are helping us in the study of the dynamics of mesoscale ocean currents. New 
products are new available, such as wavenumber spectra (Fi>, 1,83; Menard, 1983) o„ 
global maps of the statistics of the variability (Douglas et.al., 1983; Cheney 
et.al., 1983), and sample synoptic naps (Robinson et.al., 1983) of wnicj oceanogra- 
phers have to take advantage. Three further altimeter missions which are currently 
planned include the European ERS-1, the Ameri'tn TOT EX, and the French SPOT (POSEIDON 
experiment). There are many ways to evaluate the altimeter missions of these satel- 
lites, and many parameters to optimise. 

The direct measurement of sea euxfaee height is an achievement of great scien- 
tific importance and the coverage provided by a satellite data base is unique ia its 
space-time extent and along-track density. However, now and in the conceivable 
futare, sneh data will be available only during a very few costly special missions. 
Moreover, a shortcoming o/ remotely sensed data to the deep sea ocesnographe.: is the 
fact that only sarface variables are m. isared. Thus every effort mast be made to 
plan sampling strategies and analysis schemes to optimise the scientific utility of 
the obsecrations. Re are engaged ia research to sldre^s several relevant qaestioas 
along these lines. Oar approach involves simulating the acquisition and analysis of 
altimetrie data sets by performing operations on four-dimensional oceanic or 
oceanic-like fields. These fields u* obtained by either i) combining existing real 
ocean data obtained from various in situ instruments o, si) using existing observa- 
tions to drive a numerical oeet.i model (Robinson, 1984). In the first ease we fly a 
simulated altimeter over historical data and ia the second case we fly a simulated 
altimeter over simulated data. Zn oar analysis we iaclade the quantitative 
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investigation of the construction of deep ocean fields combined from both remotely 
sensed and in situ observations. Finally we mention our research on the assimilation 
of such observed fields in dynamical ocean models (Robinson £ Le si ie * 1984; Robinson 
et.al., 1984) for initialization and verification purposes and for dynamical interpo- 
lations and optimal estimates. 

Here we describe preliminary experiments; firstly of observational network 
evaluation, i.e. optimization of sampling, in order to improve mesoscale mapping and 
statistics (Section 3), and secondly the construction of cuxrent maps (Sections 2 and 

3) for the assimilation of the sea-surface height information in our quasigeostrophic 
model (Section 4) . 

2. SATELLITE ALTIMETER SIMUIATIONS IN THE POLYMOPE AREA 

In oi "*r to s ^ * Sea Surface Height (SSH) data as a satellite would do, an 
orbital model is needed. The requirements on this model are not very severe, since 
we only wart to infer statistical conclusions on the coverage, and get synoptic mid- 
ocean maps of SSH. A simple circular orbit model has been used, along with accurate 
orbital parameters of future missions (ERS-1, TOPEX, SPOT). No atmospheric drag is 
involved. The repeat cycle, cross-track distance, inclination angle, and subcycles 
are a^quately modelled. The numerical clock is accurate to 1 sec for more than 1 
year, in order to ensure exact repetitivity during long simulations. The altitude 
calculated by the model is accurate to lfc. 

Linear optimal estimation is used to restore SSH maps. This is an 'intelligent* 
interpolation method, whose estimator is statistically data-adaptive • It allows the 
integration of other remote-sensed and in situ measurements, and is consistent with 
the data assimilation techniques currently in use at Harvard (e.g. Tu, 1981; Robin- 
son, 1983K The estimation method is a space-time extension of the Gandin (1963) 
Objective Analysis method. Phase propagation is built into the signal correlation 
model. A limited set of influential points is selected at each interpolation point 
(Carter, 1983). An improved algorithm eliminates data points which are more corre- 
lated to each other than to the interpolation point. The need for such a stringent 
selection arises from the one-dimensional nature of the data. 

Restored SSH maps can be used for operational system assessment, which is 
inferred from error estimates provided by Objective Analysis (see Section 3) or by 
initializing and updating the Harvard quasi-geos trophic barociinic model (see Section 

4) . 


We have used two different data sets to represent sea-truth. The first is a 
1000 km simulation generated by the Harvard model using regionally tuned statistics 
in the POLYMODE-SDE region. The second sea-truth we have chosen is the 100 dbar 
level of the POLYMODE-SDE Mark II data set produced at Harvard. An interpolation 
scheme, similar to the one described above, has been used to combine POLYMODE-SDE 
Soviet current-meter moorings (1400-700 dbar) and XBTs (400-100 dbar), from Julian 
day 3345 (July 21, 1977) through Julian day 3725 (August 5, 1978). A mean salinity 
profile has been nsed to compute shallow geostrophic velocities relative to the 700 
dbar absolute velocity level. The grid has 281.25 km long sides and is centered on 
29.00 N, 70.00 W. 

SSH is sampled from the sea-truth along area given by the orbital model. The 
stream function in the shallow level is interpreted in terms of pressure, which is 
converted into SSH by f 0 /g. An overall mean is subtracted from the SSH time series 
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in the domain. Individual passes are assumed to be perfectly unbiased and detrended. 
As an illustration of the SSH variability in the whole POLYMODE-SDE Mark II data set, 
the following statistics are derived: 

Mean elevation = 1.41 m ( ref. to 700 dbar) 

Variance = 860 cm* 

Standard deviation = 29.33 cm 
Minimum (rel. to mean) = -28 cm 
Maximum (rel. to mean) = +33 cm 
Max. amplitude = 61 cm 

Data count (IBTs) = 1336 

The variance above includes SSH variability at time scales much larger than the 
typical mesoscale time. Thus, the standard deviation order of magnitude of 10 cm has 
to be kept in mind. 

A gaussian noise of 3 cm standard deviation is added to SSH time series. 

The SSH correlations in the Mark II data set exhibit the following scales: 

Space: SW-NE direction: 110 km Time: e-folding: 10-20 days 

SE-NW direction: 160 km 


The space correlation scale is defined as the zero-crossing of the covariances. 


In Section 3, an isotropic stationary signal correlation model has been used. 
It has the following analytical expression (see Fig. 1): 


with: 


C(r,t) 


(1 + &r -|a 3 r 3 ) 


- (ar + t l / P c t) 

e 


(l) 


a = 2.1038/R CX 
R cx = correlation scale 
Ret = correlation tine 
r = ((x-c x t) a + (y~c y t) a ) (/1 . 

Using this expression leads to a -2 power law for the kinetic energy spectroa. 
The following values have been used: 

R cx = ISO ka 

R ct - 13 days for 1000 ka rnns (error only) 

20 days for 281.25 ka rnns (with sea truth) 



Residual SSH pover spectra, as obtained by Menard (1983) or Fu (1983), along 
repeat tracks of SEASAT altiaeter data, can be used to calculate the signal correla- 
tion. If the along-track wavenuaber spectrum is denoted by E t (k ( ), the one- 
diaensional isotropic correlation function is written: 

, 4-09 

C(r) “I* | (k, ) exp (2irik,r)dk, (2) 

J -n 

i.e. the Fourier transfora of the speetrua. A typical correlation funotion for the 
aid-latitude sou:h Atlantic is given on Fig. 2. For wavelengths shorter than 300 ka. 
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the shape of this function shows good agreement with the shape of (1). 

Clearly, the Objective Analysis technique cannot be more accurate than the 
statistics used to build the estimator. The question of confidence levels on vari- 
ance and covariances is thus of critical importance. The availability of large data 
sets of SSH sampled on a homogeneous mesh on iarge space and time scales is a 
decisive advantage of satellites compared to in situ sensors: the statistics derived 
from remotely sensed data are likely to be much more reliable. In the POLYMODE-SDE 
area, estimates based upon GEOS 3 data (Douglas et al., 1983) and SEASAT data (Cheney 
et al., 1983) are available. The 3.5 year calculation of SSH variability from GEOS 3 
(Fig. 3) yields a typical value of 10 cm. The 24 day calculation from SEASAT (Fig. 

4) gives 5 cm. The discrepancy is due to considerable attenuation of the mesoscale 
signal (50-150 days) in the second case. The first value of 10 cm is more reliable, 
and in good agreement with the calculations from in situ measurements discussed 
above . 


The noise level, denoted by E, has been set to 20% of the SSH variance. The 
residual noise in the altimeter data, using the repeat track method, is expected to 
be 5 cm at the mesoscale, or less, for SEASAT and future missions: 

- Residual tropospheric error: 1-2 cm/50-500 km 

(Rain and water vapor) 

- Atmospheric loading: 3 cm/200-1000 km 

(Inverse barometer effects) 

- Residual geoid: a few oentimeters/all scales 

(Tracking errors) 

Assuming the 10 cm variability of Douglas et al., we find that E = 25%, which 
is likely to be pessimistic if the data have been correctly undersampled along 
passes. 

3. TWO EXAMPLES: ERS-1 and SPOT 


Figs 5 and 6 are RMS expected error maps for two 1000 km simulations carried out 
with regionally tuned statistics in the POLYMODE-SDE area, reap, with ERS-I and SPOT. 
The following orbital parameters are used: 

ERS-1: 14 + '/$ rotation/ day 

(Fig. 7) repeat cycle: 3 days 

cross-track distance: 910 km 

inclination angle: 98.52* 


SPOT: 14 + 5 /a< rotations/day 

(Fig. 8) repeat cycle: 26 days 

cross-track distance: 106 km 

inclination angle: 98.72* 


These two satellites stand for the two possible extremes as far as regional 
Z" . mid-latitude coverage is concerned. T0PEX (Fig. 9) is just in between (repeat cycle: 
** 10 days). ERS-1 and SPOT are sun-synchronous : the angle between ascending and des- 

)| cending passes (respectively 24.85° and 2S.33°at the equator) is sharper than for 
T0PEX (44.61° at the equator). 

.►» 

^ Let us define the overbar as the apace and time average on the maps on an 
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integer number of repeat cycles (6 for ERS-1 = 18 days, 1 for SPOT = 26 days), c r* 
(x,y,t) as the local expected error in terms of signal variance, v as the signal 
variance, h (x,y,t) and h (x,y,t) as the actual and restored SSH. The following 
statistics are derived for a 1000 km domain: 


jy IjAO A 

Global RMS expected error = * : 71.90% 
Global MS expected error = <r z 51.70% 
Global average expected error (bias) 0 


SPOT 

56.15% of st. dev. 
31.50% of var. 

0 


ERS-1 exhibits an excellent localized resolution, and a high variance of error, 
due to the absence of a subcycle to scan between adjacent tracks. The inhomogeneity 
in resolution would probably result in iigtroper model assimilation. 


SPOT shows a good general resolution, and a low error variance. The minimum 
error (29%) is higher, but the global statistics are much better. The whole error 
field propagates with the subcycle, i.e. at 25.1 cm/s, westward. Scanning between 
two successive tracks, with a subcycle of 5 days, is responsible for the good 
behavior of SPOT. TOP EX has a subcycle of 3 days, scanning eastward at 122.2 cm/s. 
SPOT is probably better for longer periods and relatively slow-moving events, while 
TOP EX is likely to perform better for energetic, rapidly evolving events st the 
sdvective time scale. Besides, TOP EX is not sunsynchronous, thus releasing the con- 
straint of an almost polar orbit. The angle between ascending and descending passes 
is more open, and the chances tc catch a mesoscale event are higher. 


Fig. 10 shows a 5-day sequence of SSH simulation and reconstruction, in the 
281.25 km domain sampled by SPOT. Also shown are the actual error field and the 
expected error field, both expressed in signal variance units. The following overall 
statistics are derived from the 26-day siamlation: 



Nark-II domain 

Global RMS expected error •= 

-r'/x 

<T Z 

ERS-1 

43.11% 

of 

s t . dev . 


Global 

MS expected error = 

T x 

18.58% 

of 

variance 


Global 

average expected error 

<k4) iA /v* , 

0 




Global 

RMS actual error » 

40.18% 

of 

st. dev. 

V 

Global 

MS actual error = 

OT^JVv : 

16.14% 

of 

variance 


Global 

average actual error = 

h - h / v'£ : 

4.95% 

of 

st. dev. 


3 


* 


The bias is clearly negligible. It la striking to see how close the expected 
and actual values are. Objective Analysis csn thus be expected to be a good method 
to interpolate along-track data, and the altimeter to catch accurately the variance, 
even in this small domain. 

The introduction of phase speeds in the correlations, and the use of a larger 
domain are both likely to lead to more definitive conclusions on the different sam- 
pling schemes. 
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4. ASSIMILATING SEA- SURFACE HEIGHTS IN QG MODELS 


Considerable recent progress has been aade in assimilating data into the Harvard 
QG open ocean model (Robinson and Leslie, 1984; Robinson et.al., 1984). Both baro- 
tropic and baroclinic experiments have been carried out in order to resolve questions 
as to the sonrce(s) of accuracy and error in Voth the objective analysis and the 
dynamical forecast (for details of the model, see: Tu, 1981; Miller et.al., 1983). 
Figs, lla-c show a barotropic experiment involving the assimilation of sea surface 
heights in the QG model. Fig. 11a is the model generated 'Veritas' data set at 
period 3.0. The objective analysis is performed on the 'Veritas' data sampled with 
the MODE-1 observational network (Fig. lib) and with the MODE-1 network and a simu- 
lated satellite track 'cross' (Fig. 11c). The addition of the satellite track 
restores much of the missing structure of the 'Veritas' data. These SSH maps are 
used to initialize the model, which computes the evolution of streamfunction using 
the discretized QG equations. There is good agreement between sea-trnth and the 
model outputs, and the westward propagation is restored. 

Fig. lid illustrates the results of a 17 day baroclinic forecast using real 
data. Note that forecast in this 144 km domain experiment maintains (with the excep- 
tion of some spikes) s level of rms difference between forecast and data of less than 
25%. The spikes have been attributed to 'unrepresentative' hydrocasts in the data 
set. This forecast experiment does not involve use of remotely sensed data or a 
satellite siaralation. This remains as a next step for the POLYMODE- SDE region. The 
Harvard model has been initialized with Rossby waves, simulated data, simulated 
satellite observed sea-surface heights (barotropic), and real data and run over flat 
bottoms and real topography. 'This model is proving to be an efficient and accurate 
component in the prediction and description of fields in open ocean regions of vari- 
ous internal dynamics' (Robinson and Leslie, 1984). 

In the suauaer of 1983 a significant step in data assimilation was accomplished 
by carrying out a real time forecast in the regime of turbulent jets and eddies in 
the California Current (Robinson et al. 1984). The dynamical model successfully 
predicted the appearance of a zonal jet in the center of the experimental region and 
also provided the means of identifying a major eddy-eddy interaction event. The 
implications of these results for the effective exploitation of satellite altimetric 
data for practical forecasts and scientific studies are substantial. 

5. CONCLUSION 

The analyses and simulations presented here show the ability of objective 
analysis to map along-traok altimeter data in a suitable way for dynamioal model ini- 
tialization, and to evaluate sampling strategies. A set of global statistical param- 
eters have been defined for that purpose. The two case studies on ERS-1 and SPOT 
show the great interest of the suboyclo. which, as a matter of fact, doubles the 
capabilities of the instrument, by being tuned on two time scales (5 days and 26 
days). Furthermore, it can be shown (Fig. 12) that the most economical way to set 
two tracks apart is to equal their distanoe to the horizontal correlation scale. This 
is also true in time with the repeat cycle and the correlation time. 5 days and 26 
days stand approximately for the adveotive and linear time scales at mid-latitudes 
respectively. The SPOT observational strategy is thus regionally adapted, as far as 
sampling parameters are concerned. It is still unclear whether in situ data, along 
with altimeter data, are eompulsory for assimilation. Different techniques of extra- 
polation on the vertical are currently tried in Harvard, in order to initialise and 
update properly the dynamical model for multi-level runs. The efficient use of these 



techniques and of those described above can be considered an essential step towards 
our knowledge of the oceanic nesoscale on large space and tine scales. 
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FIGURE CAPTIONS 


1. Correlation function for the objective analysis of altimeter data, as 
a function of distance. As an example, a 50km correlation radius has 
been chosen to draw this plot. 

2. Average correlation function from a group of repeat tracks of SF.ASAT 
altimeter data in the SV Atlantic (from Fu, 1983). 

3. Sea surface height variability in centimeters from GEOS-3 altimeter repeat 
pairs (from Douglas et al., 1983). 

4. Sea surface height variability (contour interval = 1 cm) from repeat 
tracks of SEASAT altimeter data (from Cheney et al., 1983). 

5. RMS expected error maps for a simulation of ERS-1 in a 1000 km sqnare 
domain center on 2 9N.70W (POLYMODE Soviet current meter array). Maps 
are shown for days 1,2, 3, 4. Contonr interval is 0.08. The parameter! 
of the analysis are given in the text. 

6. RMS expected error maps for a simulation of SPOT in a 1000 km square domain 
centered on 29N, 70VJ Maps are shown for days 0,5,10,15. Contour 
interval is 0.07. Note westward propagation. The parameters of 

the analysis are given in the text. 

7. Global coverage of ERS-1 (from Tavernier, 1983). 

8. Global coverage of SPOT (from Tavernier, 1983). 

9. Global coverage of T0PEX (from Tavernier, 1983). 

10. Sea-truth, restored sea-surface height, RMS actual error field and RMS 
expected error field for a simulation of SPOT in a 281.25 km square 
domain centered on 29N, 70V (Harvard POLYMODE SDE Mark II data set). 

Contour intervals are: 2 cm, 2 cm, 1% of st. dev. of SSH, 5% of 

at. dev. of SSH. The westward propagation of the southern eddy and the 
stationarity of the northern one are reproduced. 

11. Sea-truth, restored sea-truth (without satellite track), reatored sea-truth 
(with satellite traok), and results of 17 day 144 km baroclinio forecast 
using real data. 

12. IMS expected error vs. cross-track distance for an objective analysis of 
2 parallel tracks. The correlation radius was 50.0 km. Note the 
first minimum. The upward slope after the second minimum is an artifact 
of the integration in a finite domain. 
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